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FITTED AND UNFITTED DOMAIN DECOMPOSITION USING 
PENALTY FREE NITSCHE METHOD EOR THE POISSON 
PROBLEM WITH DISCONTINUOUS MATERIAL PARAMETERS 

THOMAS BOIVEAU* 


Abstract. In this paper, we study the stability of the non symmetric version of the Nitsche’s 
method without penalty for domain decomposition. The Poisson problem is considered as a model 
problem. The computational domain is divided into two subdomain that can have different material 
parameters. In the first half of the paper we are interested in nonconforming domain decomposition, 
each subdomain is meshed independently of each other. In the second half, we study unfitted 
domain decomposition, the computational domain has only one mesh and we allow the interface to 
cut elements of the mesh. The fictitious domain method is used to handle this specificity. We prove 
//^-convergence and L^-convergence of the error in both cases. Some numerical results are provided 
to corroborate the theoretical study. 


Key words. Nitsche’s method; interface problem; Poisson problem; nonconforming domain 
decomposition, unfitted domain decomposition, fictitious domain. 


1. Introduction. The two main methods that can be used for the weak en¬ 
forcement of the boundary and/or interface conditions are the Lagrange multipliers 
method and the Nitsche’s method that is a penalty based method. The Nitsche’s 
method that has been introduced in 1971 m is known to have a symmetric and a 
nonsymmetric version mill]. In this work we consider a nonsymmetric penalty free 
Nitsche’s method 012, this method can be seen as a Lagrange multiplier method, 
where the Lagrange multiplier has been replaced by the boundary fluxes of the dis¬ 
crete elliptic operator. The method does not have any additional degrees of freedom 
nor penalty parameter. 

The Nitsche’s method has been applied to nonconforming domain decomposition 
with its symmetric and nonsymmetric version by Becker et al. [3] for the Poisson 
problem. The method has been extended by Burman and Zunino m using a weighted 
average of the fluxes at the interface for the advection-diffusion-reaction problem. 
Several difficulties in the analysis can be handled by taking the right choice of weights 
(see mi 1 uni)- For unfitted domain decomposition miia the interface can cut 
elements of the mesh, we handle this problem using the fictitious domain approach 

In the second section of this paper, we study the nonconforming domain de¬ 
composition where each subdomain are meshed independently. In the third section, 
we extend the results to unfitted domain decomposition using the fictitious domain 
method, the fourth section shows a few numerical examples to illustrate the theoretical 
study. 

Let and fl 2 be two convex bounded domain in with polygonal boundary, 
these two domains share an interface P = fli n 172. We define the domain 17 = Sli u 172 
with boundary df7, an example of 17 is represented in Figure [TJ The Poisson problem 
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Fig. 1. Example of computational domain fi. 
considered can be expressed as 

fti^u fi in 12 ^ , 1 Ij 2 , 

M* = 0 on dfl n , r = 1, 2, 
u^^u^ = 0 onT, (1.1) 

/ii• ni + /i 2 Vu^ • n 2 = 0 on F, 

where rt* and fii are respectively the unknown and the diffusivity in fli, fi G (Hi) 
is a given body force. In this paper C will be used as a generic positive constant 
that may change at each occurrence, we will use the notation a < b for a ^ Cb. For 
simplicity we will write the L^-norm on a domain 0 , ||•|| 2 , 2 (e) as ||•|| 0 ■ 

2. Fitted domain decomposition. 

2.1. Preliminaries. The set (7))}^ defines the family of quasi-uniform and 
shape regular triangulations fitted to We define the shape regularity as the exis¬ 
tence of a constant Cp G R+ for the family of triangulations such that, with pK the 
radius of the largest circle inscribed in an element K, there holds 

— ^Cp ^KeTh, 1 = 1,2. 

PK 

In a generic sense we define K as a triangle in a triangulation 7)) and ■= diam(iA) 
is the diameter of K. Then we define hi := vciax.p^^'j-ihK as the mesh parameter 
for a given triangulation Tf. We study the domain decomposition problem with two 
subdomains that can be meshed independently, we make the assumption p, 2 hi ^ ^ 1/12 
and we set h := max (ft-i,/i 2 )- 

Let Vi := {v G (Hi) : u|an = 0} for i = 1,2, then u = {v},vf) G Vi x I^- ^k{K) 
defines the space of polynomials of degree less than or equal to k on the element K. 
On each domain 12^ we define the space of continuous piecewise polynomial functions 

:= {uheV-.UhUe^kiK) ViF G 7))} , k ^ 1, 

X every function in Vjf has two components Vh = We now 

recall two classical inequalities. 

Lemma 2.1. There exists Cp g R+ such that for all w G (K) and for all 
K G Th, the trace inequality holds 

MdK ^ C't Ikll^ -F hi ||Vw||^) . 
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Fig. 2. Example of mesh of fl. 


Lemma 2.2. There exists Cj e R+ such that for all Uh £ f‘k{K) and for all 
K B Th, the inverse inequality holds 

At the interface we use the notations 


M — up', 


for the jump, and 

{w} = OJiW^ + UJ2Up 
with the following weights 


{w) = UJ 2 UU^ + LOiVp 


Wi = 


UJ2 = 


h2hi 


hifl2 + ^2^1 hlh2 + h2/il 

At the interface L we define the normal vector n := ni = —n 2 , then we define 
{fiVw ■ n} = LUiftiVwi • n + uj 2 fi 2 ^W 2 • ‘n^ 


we note that 


{/rVu • n} = • ni = —• n 2 . 

To simplify the notations in the analysis we set 

_ MiM2 
hih2 + h2hi 

We now introduce a structure of patches that will be used in the upcoming inf-sup 
analysis similarly as [ail]. Let the interface elements be the triangles with either a 
face or a vertex on the interface L. We regroup the interface elements of fli in closed 
disjoint patches Pj with boundary dP^, j = Np defines the total number 

of patches. Note that P* := Let := dPj n L. For each Pj there exists two 

positive constants ci, C 2 such that for all j 
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with hi the meshsize of the subdomain considered. Let us focus on the patches 
attached to the domain fii. Each patch is associated with a function (j)j G 
defined such that for each node Xi G of the face Pj let 


(j)j {Xi) 


0 for Xi G ^i\Fj 
1 for Xi G Fj , 


( 2 . 1 ) 


with i = 1,, Nn- Nn is the number of node in the triangulation and Fj defines 
the interior of the face Fj. We define the function Vp G Vi^ such that 


Np 

vj = uj, 

1=1 

with 

Vj = Vj(j)j , Vj G R. 


( 2 . 2 ) 


(2.3) 



Fig. 3. Example Patch Pj on the interface elements, the function <f>j is equal to 0 in the 
nonfilled nodes, 1 in the filled nodes. 


In order to dehne the properties of uj, we define the Po“Projection of a function 
w on an interval I 

:= meas (/) ^ 

We may now define on each face Fj the function uj such that 

meas [Pj^ ^ f Vuj • n ds = ^ . (2.4) 

Applying the Poincare inequality, on each patch Pj the function vj has the following 
property 

3 3 


(2.5) 
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It is straightforward to observe that on each face it holds 


^ h^ ||V<||^. . (2.6) 

3 

The Lemma 4.1 of [7] allows us to write the following inequality for each patch Pj of 
the triangulation T^, 





hi "luhl 


(2.7) 


Lemma 2.3. Considering the patches Pj as defined above ^Uh e the following 
inequality holds 


Nt, 

P 

E 

j=i 




Nt 


Nt 


Nf 




F} 


Eh" pi - C'wi 2 I S/ul - CuJ2 2 


i=i 


i=i 


j=i 


M2 Vu) 


P? 


Proof. Considering the triangle inequality and the definition of the jump we can 
write 





1 ^ 


1 / 

/A 

7" HI 

1 - - -F^ 

+ 

7 I 

1 / 


7= HI 

+ 

FI 

75 


7" lu/tl - [will 


,1 .,2 72^ 


Taking the sum over the whole interface and using the triangle inequality once again 
followed by inequality (12.611 , trace inequality of Lemma 12.11 and inverse inequality of 
Lemma 12.21 we obtain 


Nt 


1 = 1 


ihuhj 

1 

/A 


p 

+ 7^E 


- 

1 j 
-< 

Nl 

+ 7^E 

Tp 2 

2 2 
Uh-K 


3 j=i 


4-1 

3 0/ —J- 



F! 0 = 1 





P 

^E 

1=1 

K 

<E 

1=1 


7" HI 


Nt 


NZ 


-f 7^1 E + 7= ^2 XI 


7" HI 


J=1 

K 

+ Cwi X 
j = l 


pil Vu 


h II p2 

1=1 " 

K 

^ -I- CuJ2 2j h2 

i j=l 


p? 


.□ 


We end this section by defining the triple norm of a function w 

2 


-E 


1 

2 

1 


+ 

fii 

7"M 


( 2 . 8 ) 


2.2. Finite element formulation. Classically for the Poisson problem (HU 
for each subdomain domain fli we obtain by integration by parts 

(y^iVu*, = (/,fo)j.j_ , Vfo G Fi. 
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By taking the sum of the boundary terms we obtain 

2 r 

V • ni,v") = lifiVu ■ n) w| ds = ({fiVu ■ n}, |z;]>p + <|mVu • n] , (n)>p , 

(HI) tells us that {fiVu • n] = 0, then we get 

2 2 
2 Vu*) - <{^Vm • n} , |n]|>p = ^ (/„ . 

i=l ' i=l 

Adding the Nitsche term, it leads to the following finite element formulation : Find 
Uh e Vj!^ such that 

Ahiuh,Vh) = LUvh) yvhSVt ( 2 . 9 ) 


where 

2 

Ah {uh, Vh) = ^ (AiiV<, - ({AiVw/i • n} , lvh}\ + • n} , [w/iDp , 

i=l 

Lh {vh) = 2 

2.3. Inf-sup stability. This section leads to the inf-sup stability of the penalty 
free scheme previously introduced, we first prove an auxiliary Lemma. 

Lemma 2.4. For Uh,Vh g with Vh = Uh + Up, Up defined by equations dza 
and (12.311 . there exists a positive constant (3 q such that the following inequality holds 

fio |||uft|||^ ^ Ah{uh,Vh)- 

Proof. Using the definition of Up, we can write the following 

K 

Ah {uh, Vh) = Ah {uh, Uh) Ah {uh, v]) . 

j=i 


Clearly we have 


and 


Ah{uh,uh)= pt^Vu 


Ml 


Ah {uh,v}) = (/riVu),,Vuj) 1 -({fiWuh ■n},v]) ^ +uJi(piWv] • n, . 

3 3 3 

Using Cauchy-Schwarz inequality and inequality dZZl) 


{fil^Uh 


,aVu])j 


> — 


pil Vu 


afj.{ 


Vu, 


j II pi 


> 


—e 
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Using the trace and inverse inequalities, (1^ and (IQ) we can write 

• n} , • n + UJ2^J^2^u\ ■ n, 

,1 


= Vul ■ n + {oJ2^J■2h2)'^ ^u'h ■ 71,072^ 

^ {{uJi^jLihiY IVuj, • n||^i -I- (ui2H2h2y ||Vu^ • a 

< {{iOi^jLihiY ||Vujj • n\\pi + {u2fJ.2h2)^ ||Vu^ • n||^i^ a 


757 ;,- 


7" luftl 



Taking the sum over the full boundary T and using trace and inverse inequalities once 
again we obtain 


K 

, aVjy^-^ 

1 

j = i 




2 e 


i=i 


7= M 




+ e {ujinihi) E II Vujj • n||^i 


7=1 




+e (W2M2/12) E ■ ”IIf 2 

J=1 " 


<—s 

i=i 


7" M 


2e 

-l-ea;2 E I ^^2 




+ ewi E 

j=i 


7=1 




Using the property (12.41) of uj we can write for each face Fj 
auji (^^lVvy ■ n, ^ a 

Using the trace inequality and inequality (I2.7p we get 


7" M 


-f auji { ^iiVv] ■ n, |uft] - |m/i]| 


awi ( AtiVuj • n, |u/j| - 


^ owi^i II • n|| 


[m/i] - luhl 


< auii^ih^ 


II II pi 




ul - { 


- Ut 


^ O- 

^ Oi 



I h, ) 

( Wi/xA = 

I h, ) 




+ 
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Taking the sum over the whole interface T and using Young’s inequality and (12.61) 




^ acoi / fiiVv] ■ n, |u/i] - fuk} 

1=1 \ / F} 

P 

E 


< 


Ca^ 

4 IT 

2 

UJlfli 

_pi 

1 1 J 

2 

CJi/ii 

2 — 

2 

2e 

7"M 

ai 

3 

Uh-K 

jpi ai 

3 

Uh-K 



Ca^ ^ 

^ — E 


< 


2e 

Ca^ ^ 


i=i 

n' 


7= Im/iI 


Wi/ii 

+ e—r— 


K 

E 

i=i 


1 — 
^h-K 


2e 


E 


7" M 




Wi/ii 

+ e~r— 

pi fti 

3 


E 

i=i 


+ ewi E Ml 2 

Fl j = i ' ^3 j = i ' 


i=i 

It allows us to write 

K K , 

^ aoji (niVvj ■ n, ^ -ewi ^ Vu^ 


2 

^h-K 

2 

F? ' 


F? 


Ni 


3 = 1 


3 = 1 


pi 


■ euj2 


E 


3 = 1 


2 

P? 


Mi 




3 = 1 


7" Iu?il 


The full bilinear form Ah now has the following lowerbound 


Mi 


Ah {uh,Vh) > Ml Vuj, + 


ni\pi 


Ml 


Mi 


+ p.+ ^i=E 

3=1 1 i=i 


n2\P' 

Mi 


+ aE 

i=i 


Mi 


Ml" 


pi 


7=K1 


pi 


with the constants 


Cq — 1 — e (2a;i + 1), 

Cb = l- 2euj2, 

Cc = a [ 1 — a— ( 5 + 


4e 


^2M1 
hi ^2 


Using Lemma 12.31 it becomes 


Mi 


A {uh,Vh) ^ ^ + (Ca - uiiCCc) 2 liMfVuj, 

Mp Ml 

+ {Cb-UJ2CC,) 2 ||m|Vu 2|P^^ + Ce 2 pM 


i=i 

2 

Mi 


2 

pi 


j=i 1 i=i 

First we fix e = min [^ 2 ( 2 Ji+i) ’ 4 i^] ■ constant Cc will be positive for 

dewi 

C (4wi + 1) ■ 


a < 
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The terms (C^ — uJiCCc) and {Ch — 0 J 2 CCc) will be both positive for 

O (e) = min \0 (1) , O O (a) = min [O (ewi), O (1)] and O (/3o) = O (q!).D 

Theorem 2.5. There exists a positive constant /3 > 0 such that for all functions 
Uh G Vjf the following inequality holds 

P\\\uh\\\ < sup --. 

Proof Considering Lemma [2.41 the only thing that we need to show is 


llki^lll ^ \hh 


The triangle inequality gives 


Nt 


ik/iiii = iiiwftiii + XI lll^llll ■ 

J=1 


The triple norm (12.81) gives 


^iir = 


+ 




Recalling the inequality (lO) and 


< 


^ IIImiiIII , 

I FJ' 


it gives the appropriate upper bound 

< luF ll|w/i||| ■ 


Using the trace inequality of Lemma dT] and the inequality (12.51) 


72 m 


< plI Vv] 


< 


\Uh\\ 


Note that O [j3) = O {jSqujI ^. 0 

2.4. A priori error estimate. The proof of the stability done in the previous 
part leads to the study of the a priori error estimate in the triple norm, the following 
consistency relation characterizes the Galerkin orthogonality. 

Lemma 2.6. If u b (f^i) x (^2) is the solution of and Uh £ the 
solution of S2.f)\} the following property holds 

Ah{u-Uh,Vh) = 0 , yvheV^. 
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Proof. We observe that Ah {u,Vh) = Lh (vh) = Ah {uh,Vh), ^Vh G □ 

In order to study the a priori error estimate, we introduce an auxiliary norm 


w + 


hi 



1 1 

III hi Vw 


2 



Lemma 2.7. Let w e (Hi) x (112) + and Vh g Vh, there exists a positive 
constant M such that the bilinear form Ah (•, •) has the property 

Ah {w,Vh) ^ M ||u;||^ |||u?i|||. 


Proof. Using the trace inequality of Lemma l2.1l and the Cauchy Schwartz inequal¬ 
ity we can show 


<{^Vwn}, 


= t 2 Vw • n + (u} 2 P- 2 h 2 ) ^ Vw • n ,72 




+ wl I • n 




<{/rVuh • n} , H>r ^ lilhlVvl ■ n ^ + ujf irlhlVvl-n ^ 7^u>| 


< w 


lilVvl 


r 

+ ^2 




7" M 


Using these two upper bound it is straightforward to show that 
2 

^ {lifVw\Vvl)^_ + ({/iVw • u}, |u?,I|>p + {{pA/vh ■ n} , |u;|>p < ||u;||^ |||u,^||| . 


i = l 


O (M) = min O (ujI ) , O (wJ )]. □ 

Proposition 2.8. If u g (Hi) x {II 2 ) is the solution of and 


Uh e Vh solution of H2.9\) . then there holds 

l||u - M^lll ^ C^h'^ |M|jpfc+i 


(H) 


with Cfj, a positive constant that depends on p and the mesh geometry. 

Proof. Let denote the nodal interpolant, the approximation properties for each 
K eTh with i = 1, 2, give 

I|w - II V (u - TT^u) \j^ + h]^ \D‘^ {u - TT^u) \u\H'‘hl(K) ■ 

Using this property and the trace inequality it is straightforward to show that 


F - ^ ||w - ^ hi K \u lH»=+i(n^) + hi K \u 


Using the Galerkin orthogonality of Lemma 1^61 the Theorem l2.5l and the Lemma [2771 
we can write 
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Using this property and the triangle inequality we obtain 


■ Uh\\\ ^ \\\u ■ 


^11 fe II 

P II llx 


The result follows. 0{Cfj) = O ^ 

Lemma 2.9. Let u e (^i) ^ (^ 2 ) be the solution of 11.Ill and Uh the 

solution of i2.9i) . then there holds 


UhWn 




'\h>=+^(q.) 


with C'^ is a positive constant that depends on p, and the mesh geometry. 
Proof. Let 2 : satisfy the adjoint problem 


—— u\ 

z* = 0 

2 ^ - 2 ^ = 0 

/iiVz^ • ni + PL 2 S /• 712 = 0 


in Hi , 7 = 1,2, 
on dPl n rii , 7 = 1,2, 
on r, 
on r. 


Then we can write 

II- <lln, = (“* - <’ 

= (V {u" -ul) ■ ni,u^ - , 

it leads to 


||u - WTillo ^ 1^^ - '^hWn, + - ^hWn^ . 

2 

^ Zi “ '^h) > li*V2*) - <{AiV2 • 77 } , |77 - 777j]|>p , 

i=l 

Ahiu- Uh, z) -2 ({flVz • 77 } , |77 - 77,j]|>p . 

Using the global trace inequality ||V 2 * • 77||p < ||'Z*||^ 2 (q.) for * = lj2, we can write 


|<{/7V2 • 77 }, |u - 77?,]>r| 

^ (^(uJi/lihi)^ ||V2^ • r7||p 4- ( 072 / 72 / 72 )^ ||V2^ • 77||p^ II 75 |77 - 77;i] 
< (^(Wi/ii)5 4- (072172)") h^ \\ 4 h^Q) ll|ii - ^h\\\ ■ 

Using Lemma [121 and ( 2 * — |r = 0 for 7 = 1, 2 
Ah (77 - Uh, z) = Ah{u- Uh, z - nlz) , 

= E(V(«^-<),li.V( 2 *-^^ 0 )n. 

2=1 

+ <{liV (2 - TtIz) • 77 } , |77 - 777,|>p , 

< ((07il7i)^ 4- (072172)^) h |z|^2(o) |||77 - 777,||| . 
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Then we get 

\\u - Uh\\Q ^ Cf, ((wi/ii) = + (a;2/x2) = ^ {h + h^)h^ |u|ff*,+i(n) ||2:|li/2(n) ■ 

We conclude by applying the regularity estimate ||^||/f 2 (Q) ^ j|M —O {C'^ = 
O [Cfi ^(wi/ri) ^ + (w2/r2) ^^ 

3. Unfitted domain decomposition. 

3.1. Preliminaries. In this section the two subdomains fli and 112 are both 
meshed with only one triangulation. Let {7^}^ be a family of quasi-uniform and 
shape regular triangulation fitted to H, the shape regularity is defined in the same 
way as in the previous section. In a generic sense, F with normal np denotes a face 
of a triangle K in the triangulation 77- Tbe mesh size is defined as h := maxKeTh^K- 



Fig. 4. Example of mesh of f2. 


Figure |4] shows an example of configuration, the mesh do not fit with the interface 
r. Let 


L!f := {iF G I iA n n, ^ 0} , 

We redefine the spaces defined in the previous section. 

U := {uGifi (L!f) :n|5n =0}, 

:= Kg 0 :«^lKGPfc(7A) ViA e , 

then := x V 2 ^. We define the set of elements that intersects the boundary 

Gh := {iA G 00 7A n r ^ 0} . 


For the sake of precision we make the following assumptions regarding the mesh 77 
and the boundary F : 
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• The boundary T intersects each element boundary dK exactly twice, and each 
(open) edge at most once for K £ Gh- 

• Let r^.h be the straight line segment connecting the points of intersection 
between L and dK. We assume that is a function of length on ; in 
local coordinates 

^K,h = : 0 < ^ < ,r] = 0} 

and 

Tic = : 0 < ^ < \TK,h\ ,r] = S {^)} . 

• We assume that for all K G Gh there exists K' ^ Gh and K r\ K' ^ 0 and 
such that the measures of K and K' are comparable in the sense that there 
exists Cq > 1 such that 


meas (K) 

'' meas {K') '' 

and that the faces F such that K n F 0 and K' n F ^ 0 satisfy 
I meas corneas {F)'^ . 

• We assume that in a triangle K intersected by the interface T, the scalar 
product between the normal of the face that does not intersects T and the 
normal of the interface n keeps the same sign in K. 

We now extend the trace inequality for v G (fl) 

ll''^llAnr ~ ^ (3-1) 

The inequality (EH) has been shown in m- Let Ei be an iL^-extension on ft*, 
Ei : F[^ {rti) Fl^ (fl*) such that = w and 

~ ll'^llR'=(ni) > fc > 0. (3-2) 


Let TT* be the standard nodal interpolant, we construct the interpolation operator tt/i 
such that TThW = tt^w) such that 






We recall the interpolation estimates forO^r^s^/c-l-l, 


2 I 


llu" - < h 


s—r L,2 








s—r — -k L,2 




yK G Th, 
yFeF. 


Using the estimate (13.21) together with (13.41) and (13.51) 


IF i//'•(n*) 


El 

FeF 


IIf'-(F) 


IffUOi) ■ 


(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 
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The stability analysis for this case is similar to the fitted case treated previously, first 
we need to adapt the structure of patches to this new configuration. Let us split the 
set Gh into Np smaller disjoint set of elements Gj with j = 1,..., Np. Let be the 
set of nodes {xn} belonging to Gj. A generic node of the triangulation is designated 
by Xn, we define the sets of nodes ij and /| such that 

Ij . \^Xn £ ^Gj I ^ ^ ^Gi A , 

Ij {^n ^ IGj I ^ ^25 ^ ^Gi A , 

now we define Pj and Pj for each Gj such that 

P/ := Gj ^{KeThl ij e K} , Pf := G, u {K e % \ /J e K} . 

Each patch Pj is constructed such that /] 0 for i = 1,2. Figure [3T] shows 

an example of two patch P^ and Pj attached to a set of interface elements Gj. 
Lj := r n Gj is the part of the boundary included in the patches Pj and Pj . For all 
j and i = 1,2, the patch has the following properties 

h < meas(rj) < h and < meas(Pj) < (3.8) 

The function (j)j attached to Pj introduced in (12.11) is modified such that 


(j)j (Xi) 


0 for Xi ^ Ij 
1 for Xi e Ij, 


with i = 1,..., Nn. Nn is the number of nodes in the triangulation T/t- Similarly as 




Fig. 5. Zoom of Figure [71 For a set Gj, left : example of Pj, the function (f>j is equal to 0 in 
the nonfilled nodes, 1 in the filled nodes / right : example of Pj. 

the fitted case we define the function vj G such that 

Np 

vj = a'^vj, (3.9) 

1=1 


with 


vj = Vjipj , Vj G R. 


(3.10) 
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The function uj has the property 

nieas(rj)~^ f Vuj • n ds = \ (3-11) 

Jr, 

It is straightforward to remark that (1211) and dm still hold for this new configuration 


Uh - ul 


<h\Vul 


'illr. 


(3.12) 


\\p. 


1 < h\\Vv]\\p, . 


(3.13) 


The Lemma 1 of [5] combined with the regularity assumptions on the mesh made 
previously allows us to extend (l2Tl) to 


nlpi ^ 


h "M 

Using the trace inequality dm, Lemma 12.31 can be extended to 


(3.14) 


Np _j 

Mm 

i=i 

We now have 


N„ 


N„ 


N„ 




r,- 


XI h"M ^ -Cuji X 


1=1 


1=1 


- Cu}2 X 

1=1 


p? 

(3.15) 


LOi = 


^J-2 
Ml + M2 


UJ2 = 


Ml 


Ml + M2 


7 = 


MiM2 


/i(Mi + M 2 )' 


3.2. Finite element formulation. To handle this problem we use the fictitious 
domain method on both subdomains Ui and U 2 . We assume mi ^ M 2 - Using the 
penalty free Nitsche’s method, the finite element formulation for the problem dm is 
written as : find Uh G such that 

Ah {uh, Vh) + Jh {uh, Vh) = Lh {vh) ^Vh G V)f, (3.16) 

where 

2 

Ah {uh, Vh) = X (MiV<, - ({fiVuh ■ n} , -f ({nVvh ■ n} , [w/iDp , 

i=l 

Lh (vh) = X {fi^Ph)n,- 
1=1 


The operator Jh is the ghost penalty [6], defined such that Jh{uh,Vh) = Jh iuhj'^h) + 

Jl K,^h) with 




Jl{ul,vl)=-/g X ^lL>lpulj,lDl^vlj)p. 
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This penalization ensures the stability in the case of small cut elements that might 
cause a dramatic growth of the condition number. The sets Tq for i = 1, 2 are often 
called the interior facets 


:= {F G G?. I F n ^ 0} , := {F G G?. | F n ^ 0} . 

is the partial derivative of order I in the direction np- Figure |6] shows an example 
of the sets J-q and J-q for a small part of the boundary. From [5D] we have the estimate 




Fig. 6. Zoom of Figurethe dashed line is the boundary F, the bold facets are in the set Tq 
left ipQ, right Tq . 




n* 


< 


' +Jl{vl,vD< ' 


3.3. Inf-sup stability. We define the norm 






i=l 


nf 




(3.17) 


Lemma 3.1. For Uh,Vh £ with Vh = Uh + v^, Vp defined by equations (1^ 
and l|3.10|) , there exists a positive constant fio such that the following inequality is 
true 


/3o|||m?i|||* < Ah{uh,Vh) + Jh{uh,Vh) 


Proof Using (13.171) we can write 


2 

Ah{uh,uh)+Jh{uh,uh) = 'f^{p.iVu\,Vu\f)Q.i+Jh{uh,Uh) S: 

i=l 


l^l^< 


+ 




{Ah + Jh) {uh,v]) = {p,iVu\,yv])p^ -{{pAJuh-n} ,v]')p 

3 ^ 

+ wi (miV nj • n, + A {uh, Vj) . 
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Using the proof of Lemmawith (13.11) (13.141) 


{liiVu\,aS/v])+aJh {uh,v]) > - 


PlnOi 




3 llpinni 


-Jh {uh,Uh)^ aJh {v],v]) 


> -e 


Ml" Vwj, 


2 Ca^^jLi II ^||2 

-- vri,- „1 

pi 4g 3 Wp. 

2 Ca^ f /ii \ I !■=— 

, - -j— 1 + — 7"M 

P/ 4e V M 2 / I 


Using the trace and inverse inequalities, (13.131) and (13.141) we have 


<[{^iVuh ■ n} ,av]') 


Co? 


3/T, ^ 2e 




+ ewi /r;fVuo +ecLi2|^2^'“ft 
r T I I p ■ 


p? 


Using (13.lip . (I3.14p . (13.121) the trace and the inverse inequalities we obtain 

1 2 1 _ 2 


awi (^iVu - • n, lu/i]' 


^ — ewi 


Ml" 


- ea;2 


M2"Vu^ 


P? 


, ^ Ca 
+ 0 ( 1 - — 


We now have the lowerbound 

1 . 2 


N„ 


Ah{uh,Vh)lf m"Vm, 


a*\pi 


+ 


M 2 " 


n*\p2 


+ 2 


1=1 




(Ca-UJiCCc) 


Mi"Vu^ 


+ {Cb - U32CCc) 


M 2 " 


p7 


+ Cc 


7"M 


with the constants 


Ca — 1 — e (2a;i + 1), 

Cb = l- 2euj2, 

Cc = all-a^(5 + ^ 
4e V M 2 


All terms are positive for e = min 


-^1 


and a = min 


4eoJi 1 11 

C(4wi-|-1) ’ 2Cj ■ 


2(2cl'i + 1) ’ 4uj2 J 

Theorem 3.2. There exists a positive constant /3 > 0 such that for all functions 
Uh G Vjf the following inequality holds 


/3|||w/illU^ sup 

Vh€V^ 


Ah iuh,Vh) + Jh {uh,Vh) 


Proof. Same as proof of Theorem 12.51 □ 

3.4. A priori error estimate. We have the following consistency relation 
Lemma 3.3. If u e (Hi) x (H 2 ) is the solution of il.ll) and Uh £ the 
solution of i3.16\) the following property holds 

Ah{u - Uh,Vh) - Jh{uh,Vh) = 0 , Mvh^V^. 
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Let us introduce the norm 


NIU = WMh + 


i. 1 — i 9 

fiih^Ww -n + -n 


Lemma 3.4. Let w g (Hi) x (H2) + and Vh £ , there exists a positive 

constant M such that the bilinear form (•, •) has the property 


Ah {w,Vh) < M ||-u;||^ |||z;?i|| 


Proof. Same as Lemma 12.71 with hi = h 2 = h. U 

Lemma 3.5. Let u g (f^i) x (^2) ctrid nh the interpolation operator defined 
by 113). For all Vh £ the following inequality holds 

Jh {T^hU,Vh) ^ Cjh'^ l|w|lH^+i(n) \\M\^ , 

withO{Cj) = 0(1). 

Proof. The continuity of u gives Jh (Ku,Vh) = 0. Using the definition of the 
interpolant nh (13.31) . the inverse inequality (12.21) . the interpolation estimate (13.6p and 
the continuity of E 

Jh {TrhU,Vh) = Jh {Tr^Eu-Eu,Vh) < h^ < h^ l|u||jj;.+i(n) ■ 

□ 

Proposition 3.6. If u e (Hi) x (H2) is the solution of il . ll ) and 

Uh £ V)f the solution of iJ.16\) . then there holds 

|||m — w/tlll ^ Cf^h Inl/i-fe+qn) ■ 

with Cffi is a positive constant that depends on p, and the mesh geometry. 

Proof. Using the orthogonality relation of Lemma 13.31 we can write 


Ah {Uh - T^hU,Vh) + JhiUh - TThUjVh) = Ah{u- TThU,Vh) - Jh{T^hU,Vh), 
applying this property with Theorem 13.21 and Lemma 13.41 we get 


0111 III ^ Ah{u-TrhU,Vh) - Jh{'!^hU,Vh) ^ .j-,, II Jh{TrhU,Vh) 

P \\\uh - T^hUlW^ < -^ M ||u - 7r,,M||^ - 


\\Vh\\ 


Fiill 


Using dSH) we can write 


\\\u - TThu\\\ < \\u - 7Thu\\^ < h \u +Ai|^ \u 

It is straightforward to observe that for any Vh £ we have IHr’iilll ^ llkiilll*, using 
this result and the triangle inequality we have 


||R - Uh\\\ < |||u - TThU\\\ + \\\uh - TThu\\\^ 

s£ |||u - TThulW + i ||m - 7r,iU||^ - 


Jh{TThU,Vh) 

llk/«IIL 


||u - M/illl ^ 1 + 


M + 1 
1~ 


h 111111/1-^+1(0) 


Applying Lemma [33] 
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Lemma 3.7. Let u e (^i) ^ (^ 2 ) be the solution of and Uh the 

solution of i3.16]) . then 

\W ~ ’^hWd ^ ^ 2 + 

with C'j^ is a positive constant that depends on fi and the mesh geometry. 

Proof. Same proof as Lemma 12.91 considering the new unfitted framework. □ 

4. Numerical verifications. In this section we verify numerically the conver¬ 
gences proven theoretically. For each case studied the domain considered is the unit 
square separated in two subdomains as it is shown in Figure [T] We use a manufac¬ 
tured solution in order to test the precision and determine the slopes of convergence. 
The manufactured solution that has been considered in this case 

u = exp(a;j/)sin(7ra;)sin(7ry). 

For fitted and unfitted domain decomposition we consider /ri = 1 and we test a range 
of value for fi 2 - 

4.1. Fitted domain decomposition. The package FreeFem-|--|- m is used to 
implement this case. We choose different values of the ratio / 11//12 and observe the 
and TJ^-error for each configuration. 




Fig. 7. Piecewise affine approximation, pi = 1, left ^ = 1, right ^ = f • 




Fig. 8. Piecewise affine approximation, pi = 1, left ^ right ^ 
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Figure [7] and [5] shows a convergence of order 0{h?) for the L^-error, this is a 
super convergence of order 0 {h 2 ) compared to the theoretical result. This super 
convergence has been observed for linear elasticity with the penalty free Nitsche’s 
method in [4]. Comparing all four graphs we observe that as the ratio /ii/ft .2 becomes 
smaller the constant becomes slightly bigger when ^2 grows. 




Fig. 9. Piecewise affine approximation, pi = 1, left ^ right ^ 




Fig. 10. Piecewise affine approximation, pi = 1, left right ^ = §■ 


We note that in all the cases the slope of convergence that has been proven theo¬ 
retically is observed as it is shown in Figures [5] and [TUI In fact the affine approximation 
considered gives slopes of convergence of order O {h) which is what has been shown 
theoretically. For hi/h 2 = 1 the meshsize are the same on both sides of F, in this 
case the influence of fi 2 is negligible, the error remains the same for every value of /12 
considered. By considering the ratio / 11//12 smaller, the nonconformity of the meshes 
on both size of gamma gets bigger but it has a very small impact on the error for the 
three nonconforming cases considered. 


4.2. Unitted domain decomposition. The package FEniCS [TU] and the li¬ 
brary CutFEM [S] have been used for these computations. 
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Fig. 11. Piecewise affine approximation, unfitted case, •yg = 0.001. 


The same super convergence is observed as in the fitted case for the L^-error. 
The case /i 2 = 1 seems to have a constant slightly smaller than the other cases. 
The H^-eTroT shows the same convergence as shown in the theory, once again, the 
difference between each case is negligible. 
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